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Abstract 

Using a coupling argument, we establish a general weak law of large 
<**] ■ numbers for functionals of binomial point processes in d-dimensional space, 

with a limit that depends explicitly on the (possibly non-uniform) density 
of the point process. The general result is applied to the minimal spanning 
tree, the fe-nearest neighbors graph, the Voronoi graph, and the sphere of 
influence graph. Functionals of interest include total edge length with 
arbitrary weighting, number of vertices of specifed degree, and number of 
^ ' components. We also obtain weak laws for functionals of marked point 

^\ , processes, including statistics of Boolean models. 

C-- ; 1 Introduction 

O 



Establishing laws of large numbers (LLN) for functionals of random Euclidean 
point sets is of considerable interest. When the point set forms the vertex set of a 
i-G ' graph, functionals of interest include total edge length with arbitrary weighting, 

number of edges, and number of components. Relevant graphs include those 
in computational geometry, such as the minimal spanning tree, the fc-nearest 
neighbors graph, Voronoi graph, and sphere of influence graph (these graphs 
are formally defined in Section 2.3). When the random Euclidean point set 
is a marked point set, then functionals of interest include those arising in the 
stochastic geometry of Boolean models. 

For many functionals, subadditivity works well as a basic tool; see [27] and 
[29] for surveys. This is the case with power-weighted edge length functionals. 
For example, if G(n, f) is the minimal spanning tree (MST) on n i.i.d. random 
d- vectors with common density / on R d , and if |e| denotes the length of the 
edge e, then the following asymptotics hold for the sum Y^ e eG(n f) M" °f tne 
power-weighted edge lengths: 

Department of Mathematical Sciences, University of Durham, South Road, Durham DH1 
3LE, England: mathew.penrose@durham.ac.uk 

Department of Mathematics, Lehigh University, Bethlehem PA 18015, USA: 
joseph.yukich@lehigh.edu 

1 Research supported in part by NSA grant MDA904-01-1-0029 



Theorem 1.1 If either (i) 1 < a < d, J Rd f{x)^ d a ^ d dx < oo, and J Rd \x\ r f {x)dx < 
oo for some r > d/(d — a) or (ii) a > d and f has support on [0, l] d , and is 
bounded away from zero, then as n — > oo 



E 

eeG(n,/) 



n { a -d)/d Y \e\ a ^C(a,d) f{x) {d - a ^ d dx a.s., (1.1) 



where C(a, d) is a positive constant. 

For a proof see [29] for case (i) and see [30] for case (ii). 

However, many functionals are not amenable to subadditive methods. For 
example, if the edge lengths |e| in a Euclidean graph G are weighted by a gen- 
eral function (f>, giving a sum of the form X) e eG 'KM)' then subadditive methods 
break down, and the LLN behavior is much less well understood. Functionals in- 
volving the Voronoi, Dclaunay, and sphere of influence graphs are also generally 
not amenable to subadditive methods. Functionals of marked point processes 
are usually not subadditive cither. Moreover, even when subadditive methods 
are applicable, they provide little information on the numerical values of limiting 
constants such as C(a, d) appearing in (1.1). 

An alternative approach is the so-called "objective method" . Steele [27] 
coined this term for a philosophy whereby, loosely speaking, one uses the locally 
Poisson nature of a binomial point process to describe the limiting behavior of 
functionals on finite point sets in terms of related functionals defined on infinite 
Poisson point sets. Aldous and Steele [1] used this idea to analyze certain 
functionals associated with the MST on uniform points, but one might expect 
it to be applicable to any functional, including those defined over non-uniform 
point samples, consisting of contributions which are locally determined in some 
sense. As noted in [1, 27], making formal sense of this intuition is not always 
quite so simple as one might imagine. 

In an attempt to formulate in general terms the idea of locally determined 
contributions, Penrose and Yukich [23, 24] introduce a concept of "stabilizing" 
functionals and essentially use the objective method to establish a strong law 
of large numbers (Theorem 3.2 of [24]) for stabilizing functionals on uniform 
point sets in M. d . Jimenez and Yukich [15] obtain sufficient conditions yielding 
laws of large numbers for sums involving general edge weights and non-uniform 
point sets, but their conditions are rather strong and are limited to increasing 
functions </>. 

The goal of this paper is to use the objective method to provide relatively 
simple conditions guaranteeing a general weak LLN for stabilizing functionals 
on possibly non-uniform point samples of size n. We illustrate the diverse ap- 
plications of the general LLN by obtaining weak laws for functionals of spatial 
point processes in computational geometry as well as functionals of marked 
point processes, including those arising in packing processes and the stochastic 
geometry of Boolean models. 

Many stabilizing functionals arc defined in terms of graphs G which are 
themselves stabilizing, i.e., locally determined in a sense to be made precise 



below; stabilizing graphs include the MST, fc-nearest neighbor, Voronoi, and 
sphere of influence graph. Given a stabilizing graph G, the theory applies to 
functionals such as the number of leaves, the number of components, and the 
sum of weighted edge lengths J2 e eG <^(l e l)- ^ n these graphs, edges are between 
"nearby" points, and since the density of points grows in proportion to n, the 
typical distance between nearby points can be thought of as decreasing in pro- 
portion to rT x l d . Therefore we consider sums of the form X^eeG ^("^M)' 
establishing weak LLN behavior. The limiting constants are defined explicitly 
in terms of </>, the density /, and certain graphs on Poisson processes, thereby 
providing extra information on the value of limiting constants such as that aris- 
ing in (1.1) that is not given by subadditive methods alone, even in the classic 
case when <p is the identity function. 

2 Main Results 

2.1 Terminology 

In Section 2.2 we shall formulate a collection of general LLN results. Before 
doing so we need some terminology. 

Given X C M. d and a positive scalar a, let aX := {ax : x G X}. Given y G M d 
set y + X := {y + x : x G X}. For x G M. d , let |x| be its Euclidean modulus and 
for r > 0, let B(x;r) denote the Euclidean ball {y G M. d : \y — x\ < r}. Let 
denote the origin of M. d . 

Suppose £(x; X) is a measurable M + -valued function defined for all pairs 
(x, X), where X C M. d and x is an element of X. We assume that X is locally 
finite, i.e., contains only finitely many points in any bounded region. Suppose £ 
is translation invariant, i.e. £(y + x; y + X)= £_(y; X) for all y G M. d and all x, X. 
Then £ induces a translation-invariant functional H^ defined on finite point sets 
X C M. d by 

^(^):=^^;4 (2-1) 

xex 

Functionals admitting the representation (2.1) include the total edge length, the 
total number of edges, the total number of components, and total number of 
vertices of fixed degree of Euclidean graphs. Later on, in cases with x <£ X it 
will be useful to abbreviate the notation £(x; X U {x}) to £(x; X). 

We probe the behavior of the functional H^ by inserting an extra point into 
its domain. For "typical" point sets X, it is conceivable that the contribution 
£(x; X) is not affected by changes in X which are far from x. We formalize 
this notion as follows. For any locally finite point set S C M. d , and any integer 
m G N define 



and 



£(S; m) := sup(ess sup m .„{£(0; (S n B(0; m)) U A)}) 

nSN 



£(«S; m) := inf (ess inf m> „{£(0; (S n B(0; to)) U A)}), 

— nGN 



where ess sup m „ (respectively ess inf TO) „) is essential supremum (infimum), 
with respect to Lebesgue measure on R dn , over sets ylcR' i \ -6(0; m) of cardi- 
nality n. Define £00 (5), called the limit of £ on 5, by 

£00 (5) := lim sup £(5; m). 

m — >oo 

We shall say the functional £ stabilizes on 5 if 

lim |(5; m)= lim £(5; to) = £00 (5). (2.2) 

m — >oo m — >oo — 

For t G (0, 00), let V T be a homogeneous Poisson point process of intensity 
r on M. d . We are interested particularly in functional that stabilize almost 
surely on V T - Note that with probability 1, £,(V T ', m) is nonincreasing in to and 
£,( r P T ;m) is nondecreasing in m, so they both converge. Stabilization means 
they converge to the same limit, almost surely. The present formulation of 
stabilization is weaker than that of [24]. Any functional £(#; X) which depends 
only on the points of X within a fixed distance of x is stabilizing on V T . 

We are interested in functionals on spatial point processes involving non- 
uniform points, defined as follows. Let X\,Xi, ... be i.i.d. <i-dimcnsional ran- 
dom variables with common density /, which is fixed but arbitrary. Define the 
induced binomial point processes 

X n := X n (f) := {X u ...,X n }, n £ N. (2.3) 

Our general limit theory is not for H^(X n ), but for H^ n (X n ), where we define 

^ n {x;X):=^ d x;n^ d X). (2.4) 

To obtain a LLN for H^ n (X n ) we use the following approach. By coupling 
n 1 ' d X n to a Poisson process of varying intensity, we show that the local behav- 
ior of ^{n 1 ' d Xi;n}^ d X n ) is approximated by the local behavior of the coupled 
Poisson process. If the functional £ stabilizes on homogeneous Poisson point 
processes, then a conditioning argument shows distributional convergence of 
^(n^-^X^n^-^Xn) for each X, £ X n . Under appropriate moment conditions on 
£, this gives a weak LLN for H^ n . This formalizes the intuitive notion that the 
limiting behavior of H^ n on finite sets is related to the behavior of £oo on the 
infinite set V T . 

Many of the applications that we consider are concerned with functionals of 
graphs of the form G — G(X) defined for each locally finite point set X C\M. d , 
where either G(X) or (in the case of the Voronoi graph) its planar dual has 
vertex set X. 

We shall say G is translation invariant if translation by y is a graph isomor- 
phism from G(X) to G(y + X) for all y £ M. d and all locally finite point sets X. 
We shall say G is scale invariant if scalar multiplication by a induces a graph 
isomorphism from G(X) to G(aX) for all X and all a > 0. 

It is useful to have a notion of stabilization for these graphs. Given G, and 
given a vertex x £ X, let £(x; G(X)) be the set of edges of G(X) incident to x 



(or for the Voronoi graph, the set of edges whose planar duals are incident to 
x). Let V T fi := V T U {0}. We shall say that G stabilizes on V T if there exists a 
random but almost surely finite variable R such that 

£(0; G{V T ,o)) = £(0; G{V Tfi n B(0; R)) U A) 

for all finite A dR d \B(0;R). 

Stabilization of the graph G says that the local behavior of the graph in a 
bounded region is unaffected by points beyond a finite (but random) distance 
from that region. As we shall see, the minimal spanning tree and the fc-nearest 
neighbors, Voronoi, Delaunay, and sphere of influence graphs are all stabilizing 
on V T , t e (0, oo). 

2.2 General LLN results 

The following theorem places the objective method in a general context, shows 
that the asymptotic behavior of H^ n (X n ) is sensitive to the underlying density 
/, and explicitly identifies the asymptotic constants in terms of / and the limit 
functional £oo. It will be proved in Section 3. 

Theorem 2.1 (General LLN) Suppose q — 1 or q = 2. Suppose £ is almost 
surely stabilizing on V T , with limit £ooCPt), for all r € (0, oo). If £ satisfies the 
moments condition 

8avE[t(n 1 ' d X 1 ;n 1 / d X n ) p ] < oo, (2.5) 

nEN 

for some p > q, then as n — > oo, 



n 



-i 



H in (X n )^ E[i 00 {V f{x) )]f{x)dx iniX (2.6) 



Obviously, for there to be any possibility at all for the mean of the left side of 
(2.6) to converge to a finite limit, the moments condition (2.5) must hold for 
p = 1. In this sense, when q = 1, the moments condition (2.5) is close to being 
the best possible. 

A simplification arises in the case where there is a constant 7 > such that 
£ satisfies the relation 

£,(ax;aX) = a?^{x; X) 

for all positive scalars a and all finite point sets X and x G X. In this case 
we say £ is homogeneous of order 7. Homogeneity of order 7 implies that 
£„(x; X) — n 7 / d £(x; X). Moreover, almost sure stabilization on V\ with limit 
£00 (Pi), together with homogeneity of order 7, implies stabilization on V T with 
limit T~ 1 ^ d S, oa (Pi), for any r > 0. Therefore the L q limit in (2.6) simplifies to 

E[UCPi)] I f(x) (d -^ /d dx. (2.7) 



Even simpler is the special case where the function £ is scale invariant, i.e., is 
homogeneous of order 0. In this case the expresion (2.7) simplifies to -E[£ o(?'i)], 
and the large n behavior of scale invariant functionals is not sensitive to the 
density of the underlying point set. 

Theorem 2.1 admits the following extension to functionals defined on marked 
point sets. Let (/C, Jx, Pfc) be a probability space. A marked point set is a subset 
of R d x /C, to be denoted X where X is an (unmarked) subset of R d and the 
tilde indicates that each element x of X carries a mark in /C (and with its mark, 
is denoted x). In this context, a functional £(x;X) is said to be translation 
invariant if for all y £ M. d , and for any element x of any marked point set X, we 
have £(t v (x); t v (X)) — £(#; X) where t v is the translation operator sending any 
element (x,t) € R d x K, to (y + x, t) (i.e., leaving the mark unchanged). 

In the random setting, assume the marks are i.i.d. with distribution Pjc- We 
are interested in the cases where X is the point process X n or V T ', in both cases 
assume the mark values are independent of the set X. We say that £ stabilizes 
on the marked Poisson point process V T if (2.2) holds with S replaced by P T . 

If £ is translation invariant and almost surely stabilizing on V T , with limit 
£,oo(jPt)i t G (0, oo), and if £ satisfies the moments condition (2.5) for some 
p > q, then as n — > oo, we obtain a version of (2.6) for marked processes, 
namely, 

n- l H in {X n )^ f E^ 00 (V f(x) )}f(x)dx inL«. (2.8) 

In applications of (2.8), it will be clear that we are considering marked point 
processes and for simplicity we will thus suppress mention of the tilde. 

Many applications of Theorem 2.1 arc defined in terms of functionals of 
graphs arising in computational geometry. Suppose the graph G := G(X) is 
defined for all locally finite X. Given G, functionals such as total length, number 
of edges, or number of edges less than some specified length are of interest. These 
and other functionals may be interpreted as a total of (^-weighted edge lengths, 
i.e., as a sum 

eeG(X) 

with (j) : [0, oo] — > [0, oo) a specified function. Also of interest are the number 
of components, which we denote K (X), and, for any specified finite connected 
unlabeled graph T, the number of vertices x G X for which G(X) contains a 
subgraph isomorphic to T with a vertex at x, which we denote Vp (X). Let uq 
be the order of the component (i.e., the number of vertices in the component) 
containing the origin of G(Pi y o), and let Ey be the event that G{V\fi) contains 
a subgraph isomorphic to T with a vertex at the origin. 

Note that both K G (X) and V r a (X), as well as L^(X) in the case <j> = 1 (the 
total number of edges) are scale invariant functionals of X. 

The following general result, proved in Section 3, is a consequence of Theo- 
rem 2.1. 



Theorem 2.2 Suppose G is translation and scale invariant and stabilizes on 
V\ . Then 

rT 1 K G {X n ) -► E[a G x \ in L 2 (2.10) 

and for any finite connected graph V , 

n _1 Vr G (#„) -> P[Er] in L 2 . (2.11) 

Moreover, if q = 1 or q = 2 and <fi : [0, oo] — > [0, oo) is a specified function with 



supE 




< oo 



(2.12) 



for some p > q, then 



n- x L G An x l d X n ) 



* E ■ 

ee£(0;G(7>i, )) 



/W 1 ' 



/(a;)da; in L 9 . (2.13) 



The conclusions of Theorems 2.1 and 2.2 can be strengthened in more than 
one way. If, for example, q > 2 is an integer, and if £ satisfies the moments 
condition (2.5) for some p > q, then a modification of the coupling arguments 
given in Section 3 yields convergence in L q . Moreover, for many functionals, 
convergence of the means E \n~ 1 H^ n (X n )\ , as given by Theorems 2.1 and 2.2, 
implies almost sure and even complete convergence using concentration inequal- 
ities involving either isoperimetric methods or sums of martingale differences; 
we will not pursue this but refer to [15, 24, 29] for details. We will also not seek 
to establish rates of convergence in our general results. 

Remarks. 

(a) In the special case where <fi(x) = x a (i.e., power- weighted edges), the left- 
hand side of (2.13) equals n {a - d ^ d L G (X n ), while the right-hand side of (2.13) 
simplifies to C G (d,a) J Rd f(x)^- a ^ d dx, with C G (d,a) := |££ e <K|e|), the 
sum being over all e G £(0; G('Pi.o))- This provides extra information about, 
for example, the limiting constant in (1.1). 

(b) The stabilizing hypothesis of Theorem 2.1 can be weakened to one requir- 
ing that £ be almost surely stabilizing on Vf( x ), with limit £,oo{'Pf( x )), for almost 
all x in the support of the density /. These weakened hypotheses are used in 
Theorem 2.7 below. The moments condition (2.5) is not always easy to check, 
but is obviously true for any p when the functional £ is uniformly bounded. In 
Section 2.3 we shall verify the moments condition (2.12) for various graphs. 

(c) Even without the moments condition (2.5), the stabilization assumption 
in Theorem 2.1 is enough to guarantee that the right side of (2.6) is a lower 
bound for lim inf E[n^ l H^ n (X n )] ; this is proved by following the proof of Lemma 
3.2 below and applying Fatou's Lemma. Likewise, even without the moments 
condition (2.12), a weaker version of (2.13) holds in which the right hand side 



is a lower bound for the liminf of the expectation of the left hand side. With 
the moments condition (2.5) (respectively, (2.12)) the integral on the right hand 
side of (2.6) (respectively, (2.13)) is finite. 

(d) Limit laws such as (2.10) and (2.11) for scale invariant functionals are 
of interest in multi-dimensional scaling [5]. In this context, given the matrix of 
interpoint distances between pairs of points in a data set X, one seeks to identify 
the dimension in which a data set lives. To identify the underlying dimension, 
it is useful to study scale invariant functionals of the interpoint distances, since 
these are precisely the functionals whose asymptotics are sensitive only to the 
dimension of the support of the distribution of the data points and not on the 
underlying density. This approach to dimension identification may be relatively 
inexpensive from a computational point of view [6] . 

(e) Without further conditions on £ and H^ n , we are unable to obtain asymp- 
totics for H^ n over point sets consisting of random d- vectors having a law with 
a singular component. 

(f) Suppose the conditions of Theorem 2.1 are satisfied with q = 2. Since £ 2 
is almost surely stabilizing on V T , with limit ^(Vt), and since £ 2 satisfies the 
moments condition (2.5) for some p > 1, it follows that the sample variance of 
{£ n (Xf, X n ), 1 < i < n}, namely the quantity 



converges in L 1 to Var(^ oc (7 :, /(x 1 )))- 

(g) The limit (2.11) says, loosely speaking, that the number of vertices of 
G(X n ) satisfying any property determined by the local graph structure within 
a bounded graph distance exhibits LLN behavior. For example, it yields a LLN 
for the proportionate number of vertices of any fixed degree, among other things. 

(h) A version of Theorem 2.2 also holds for directed graphs. In this context 
the limit (2.13) holds without the factor of 1/2 and with £(x; G(X)) defined to 
be the set of edges going into x. 

2.3 Applications 

The applications of Theorems 2.1 and 2.2 range from the treatment of function- 
als in computational geometry to the statistics of Boolean models as well as to 
an analysis of packing processes. The following discussion is not meant to be 
an exhaustive treatment of applications, but is merely meant to indicate the 
variety of uses of the main theorems. 

(a) Minimal spanning tree. Given a locally finite set X C M. d , d > 2, 
let MST(,Y) be the graph with vertex set x obtained by including each edge 
e := {x, y} for which there is no path in X from x to y consisting of links 
which are all shorter than e. If X is finite with distinct inter-point distances, 
then MST(A') is the minimal spanning tree on X, i.e. the connected graph with 
vertex set X of minimal total edge length; see Alexander [2] . If X is infinite, then 



MST(A') is the so-called minimal spanning forest on X [2]. Clearly MST(A') is 
translation and scale invariant. 

The following theorem considerably expands upon Theorem 1.1. We say 
that 4> has polynomial growth of order a < oo if <p(x) < C(l + x a ). 

Theorem 2.3 Let G{X) = MST(A'). Then (2.11) holds for any finite con- 
nected r. Also, (2.13) holds with q — 2 if either (i) (f> is bounded, or (ii) <f> has 
polynomial growth of order a < oo, the support of f is a convex polyhedron and 
f is bounded away from infinity and zero on its support, and (2.13) holds with 
q = 1 if (in) <j> has polynomial growth of order a < d and J Rd \x\ r f(x)dx < oo, 
for some r > max{ad/(d — a), d/(d — a)}, and J„ d f(xy d ~ a >' d dx < oo. 

Special cases of Theorem 2.3 include an L 2 version of Theorem 1.1 (take <j> 
to be a power function) and the L 2 LLN for the empirical distribution function 
of rescaled edge lengths in the MST, first obtained by Penrose [21] (take <f>(x) := 
lr 0j t](a;)). For uniformly distributed points, Bczuidcnhout et al. [4] obtained an 
almost sure version of the latter result. The first part (2.11) of the conclusion 
in Theorem 2.3 yields an L 2 version of the results of [28] and Theorem 5.2.2 of 
[27] , concerning the number of vertices of fixed degree. 

Proof of Theorem 2.3. The set of edges of MST('Pi.o) incident to the origin 
is unaffected by any additions or deletions of points outside a ball of random 
but almost surely finite radius, i.e. G(X) = MST(A') is stabilizing on V\. This 
follows from the definition of MST(A') and the simultaneous uniqueness of per- 
colation clusters (Alexander [3]), and the almost sure existence of a "blocking 
set" of Poisson points lying in some annulus surrounding the origin and preclud- 
ing the possibility of any edge from the origin to the exterior of that annulus. 
For more details see Lee [16]. 

Theorem 2.2 then gives us (2.11). To prove (2.13) it remains only to verify 
the moments condition (2.12). Note first that points in the MST have a degree 
which is uniformly bounded by some finite constant C(d) (Lemma 4 of [1]). 
Hence under condition (i), i.e. if <f> is bounded, then ^2 e ^£( x -G(x)) 'KM) ^ s a ^ so 
bounded for any point set X and x € X, and thus (2.12) is clearly satisfied, 
showing that (2.13) holds in this case. 

The uniform bound on vertex degrees implies that for any p > 1, there is a 
second constant C(p, d) such that for any X and any x <E X, 

E <M™ 1/d M)j <C(p,d) J2 <f> P (n 1/d \e\). (2.14) 

K ee£(x;G(X)) J ee£(x;G(X)) 

Suppose condition (ii) holds, i.e. <f> has polynomial growth of order a < oo 
and the support of / is a convex polyhedron with the restriction of / to its 
support bounded away from zero and infinity. Then by the uniform bound on 
vertex degrees, the right hand side of (2.14) is bounded by a constant plus a 
constant multiple of 

n ap ' d max |e| ap . (2.15) 

ee£(x;G(X)) 



Lemma 2.1 of [30] shows that 

£[max{|e| ap : e g £(X i; G(X n ))}} = 0( n - ap / d ) 

and therefore the left side of (2.12) is uniformly bounded by a constant. There- 
fore, under condition (ii), the conclusion (2.13) (with q = 2) holds. 

Finally consider case (iii). Put £ n {x; X) '■= J2ee£(x G(X)) 0(« 1 ^ d |e|). If 4> nas 
polynomial growth of order a then by exchangeability and (2.14), 

n 

E[(UXi;X n ))P}=n- 1 Ej2(UX l ;X n )r 

n 

<C(p,d)n- l EY, E r(n 1/d \e\) 

»=1 ee£(Xi-,G(X n )) 

= 2C{p,d)n- l E Yj 4> P {n 1/d \e\) 
e£G(X n ) 

<C 1 (p,d) + C 2 (p,d)n^P- d ^ d E J2 M° P - ( 2 - 16 ) 

e£G(X n ) 

We now bound the right hand side of (2.16). Let L a P(X) := J2 e eG(X) \ e \ ap 
and for all k > 1, let A k denote the annular shell centered around the origin of 
M. d with inner radius 2 fc and outer radius 2 k+1 (Let Ao be the ball centered at 
the origin with radius 2). Note that, as in (7.21) of [29], 

L ap (X n )< Y L a P(X n nA k )+C(p) max \Xi\ a ", (2.17) 

z — * l<i<n 

0<k<s(n) 

where s(n) is the largest k such that Ak f~l X n is not empty We need to show, 
after taking expectations and dividing by n( d ~ ap ^ d , that the two terms on the 
right arc bounded uniformly in n. By Jensen's inequality and the growth bounds 
L a P(X) < C(diam,Y) a P(card(A'))( d - a P)/ d (see Lemma 3.3 of [29]), the first term 
is upper bounded by 



Cj22 kap (P[X 1 eA k ] 



\{d—ap)/d 
fc>l 



which is finite by the intcgrability hypothesis (see p. 85 of [21]). The second 
term is bounded by 



C( P ) / P 



max IXXP >tn^ d - ap ^ d 

Ki<n 



dt (2.18) 



< nC(p) / P hxtfrtW-od > jA/W-«p) J dt (2.19) 



o 



By Markov's inequality together with the moment condition J Rd \x\ r f (x)dx < 
oo, for some r > ad/(d — a), this last integral is finite. Thus under condition 
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(iii), we can choose p > 1 so that (2.16) is uniformly bounded, and hence the 
moments condition (2.12) is satisfied, showing the validity of the conclusion un- 
der condition (iii). | 

(b) k-nearest neighbors graphs. Let A: be a positive integer. Given a locally 
finite point set X C M. d , the /c-nearest neighbors (undirected) graph on X, 
denoted NG(Af), is the graph with vertex set X obtained by including {x, y} as 
an edge whenever y is one of the k nearest neighbors of x and/or x is one of 
the k nearest neighbors of y. The fc-nearest neighbors (directed) graph on X ', 
denoted NG'(X), is the graph with vertex set X obtained by placing an edge 
between each point and its k nearest neighbors. If the k-th nearest neighbor of x 
is not well-defined (i.e., if there is a "tie" in the ordering of intcrpoint distances 
involving x), use the lexicographic ordering as a "tie-breaker" to determine the 
k nearest neighbors. Such a tie has zero probability for the random point sets 
under consideration here. 

It is clear that the /c-nearest neighbors graphs are translation and scale 
invariant. The following generalizes and extends the asymptotics for the sum of 
power- weighted edge lengths by McGivney [18], Yukich (Theorem 8.3 of [29]), 
and Jimenez and Yukich [15], who limit attention to continuous densities and 
increasing <f>. It also extends Theorem 2 of Eppstein, Paterson and Yao [8] , who 
prove convergence of the mean number of components for uniform samples. 

Theorem 2.4 Let G(X) denote either NG(Af) or NG'(Af). Then (2.10) and 
(2.11) hold. Moreover, (2.13) holds, with q = 2, if either (i) <fr is bounded, 
or (ii) <fi has polynomial growth of order a < oo 7 the support of f is a convex 
polyhedron and f is bounded away from infinity and zero. Finally the directed 
graph version of (2.13) holds for NG (X) with q = 1 if (p has polynomial growth 
of order a < d, L d f(xy->'dx < oo, and L d \x\ r f(x)dx < oo for some 
r > d/(d — a). 

Proof. We apply Theorem 2.2. As shown in Lemma 6.1 of [23] (even though 
the definition of stabilization there is slightly different), the set of edges incident 
to the origin in NG(Pi.o) is unaffected by the addition or removal of points 
outside a ball of random but almost surely finite radius, i.e. the graph G(X) = 
NG(X) is stabilizing on V\. Similar arguments show that NG'(.Y) is stabilizing 
on Vi. Therefore Theorem 2.2 yields (2.10) and (2.11) for NG(<Y) and NG'(AT). 

To prove (2.13) it remains only to verify the moments condition (2.12). This 
condition is checked, under any of the conditions (i), (ii), or (iii), in very much 
the same manner as for the MST of Theorem 2.3. The existence of a uniform 
bound on the degree of vertices in the nearest neighbors graph is Lemma 8.4 of 
[29]. Under condition (iii), we note that since NG' is subadditive without any 
error term, the last term in (2.17) is not needed, eliminating the need for the 
condition r > ad/(d — a). | 

Hcnzc [12] considers the fraction of points in X n which are the nearest neigh- 
bors of exactly j other points; he also considers the fraction of points in X n that 
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are the Ith nearest neighbors to their own fcth nearest neighbor. He shows that 
the limits of these fractions converge to explicit but rather complicated limiting 
constants. A directed-graph version of (2.11) in Theorem 2.2 yields Henze's 
results and shows that one can interpret his limiting constants in terms of a 
functional evaluated at a point in the origin of the homogeneous Poisson point 
process. 

(c) Voronoi and Delaunay graphs. Given a locally finite set X C M. d , and 
given x € X, the locus of points closer to x than to any other point S X is called 
the Voronoi cell centered at x. The graph on vertex set X in which each pair 
of adjacent cell centers is connected by an edge is called the Delaunay graph on 
X; if d — 2, then the planar dual graph consisting of all boundaries of Voronoi 
cells is called Voronoi graph generated by X. Edges of the Voronoi graph can 
be finite or infinite. Let DEL(A') (respectively VOR(A')) denote the collection 
of edges in the Delaunay graph (respectively, the Voronoi graph) on X. The 
Voronoi and Delaunay graphs are clearly scale and translation invariant. 

Theorem 2.5 Let d = 2 and let G(X) = VOR(Af). Then (2.11) holds. Also, 
if 4> has polynomial growth with 4>(oo) = 0, and if the support of f is a convex 
polygonal region and if f is bounded away from infinity and zero on its support, 
then (2.13) holds with q = 2. 

This result adds to McGivney and Yukich [19] and Jimenez and Yukich [15], 
who require continuous / and functions 4> which are either power functions or 
increasing. If cj> is the identity and the support of / is the unit square, then the 
right hand side of (2.13) simply reduces to 2 J, Q ^ 2 f{x) 1 / 2 dx (see Theorem 1.2 
of [19]). 

Proof. Once again we apply Theorem 2.2. As in [23], we can verify G(X) 
is stabilizing, i.e., the Voronoi cell centered at the origin for "P\ t o is unaffected 
by changes beyond a random but almost surely finite distance from of the ori- 
gin. Let £ (X; VOR(Af)) denote the edges of the cell around X\ in VOR(A'„). 
Concerning the moments condition (2.12), a modification of Lemma 8.1 of [23] 
shows that under the prescribed conditions on / and its support, 



supE 



.ee£(Xi;VOR(Ar n )) 



(n 1/d \e\r 



< oo (2.20) 



for all a > and p > 2. A modification of the proof of Lemma 2.6 of [19] 
(pp. 286-87) shows that the cardinality of £(Xi;VOR(X n )) has a finite pth 
moment for all p > 2. Combining this with (2.20) thus shows that for all <fi with 
polynomial growth 



supi? 
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v ee£(Xi;VOR(A-„)) 
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< oo (2.21) 



for all p > 2. Thus (2.12) holds and the conclusion (2.13) follows. | 

A similar result to Theorem 2.5 holds if G(X) is taken to be the Dclaunay 
graph and £(x; DEL(<Y)) denotes the edges in DEL(.Y) incident to x. In this case 
the conclusion (2.11) yields a LLN for the number of vertices of the Delaunay 
graph of fixed degree to, for any to = 3, 4, 5, .... Since this quantity is the same 
as the total number of cells in the Voronoi graph on X n which are m-gons, 
this adds to results of Hayen and Quine [11] who determine the proportion of 
triangles in the Voronoi graph on V\. 

Finally, for each t > consider the case where £(x; X) equals 1 or according 
to whether the area of the Voronoi cell around x is bounded by t or not. This is 
one case where it is natural to use Theorem 2.1 rather than Theorem 2.2, and 
that result yields a LLN for H in (X n ) := ]T" =1 ^(n 1 / d X i ; n x l d X n ), i.e., a LLN for 
the empirical distribution function of the rescaled areas of the Voronoi diagram 
on X n . 

(d) Sphere of influence graph. Given a locally finite set X C R d , the sphere of 
influence graph SIG(A') is a graph with vertex set X, constructed as follows: for 
each x s X let B x be a ball around x with radius equal to Tah\ ye x\{x}{\y ~ x \}- 
Then B x is called the sphere of influence of x. Draw an edge between x and y 
iff the balls B x and B y overlap. The collection of such edges is the sphere of 
influence graph (SIG) on X . It is clearly translation and scale invariant. 

The following LLN is apparently new, even for the identity function <f>(x) = 
x. In the case 4>(x) = 1 it extends a result of Furedi [9] on the mean number of 
edges of the SIG on uniform point sets (Furedi identified the limiting constant 
in this case). 

Theorem 2.6 Let G(X) = SIG(Af). Then (2.10) and (2.11) both hold. If <f> 
has polynomial growth and if the support of f is a convex polyhedron and if f 
is bounded away from infinity and zero on its support, then (2.13) holds with 
q = 2. 

Proof. We apply Theorem 2.2. As in [23], we can check that the edges incident to 
the origin in SIG (Pi, o) are unaffected by changes beyond a random but almost 
surely finite distance of the origin, i.e. the graph G{X) = SIG(A') is stabilizing. 
Concerning the moments condition (2.12), the arguments of [23] (Theorem 7.2) 
show that for any a > 0,p > 1, we have under the prescribed conditions on /, 



sup_E 

raSN 



^ (n 1/d Mr 



.ee£(Xi;SIG(Af n )) 



< oo. (2.22) 



Moreover, since the third moment of the degree of vertices in the SIG on X n 
is uniformly bounded over all vertices (see e.g. pp. 142-43 of [13]) we obtain for 
any <j) with polynomial growth 

p-\ 

.ee£(Xi;SIG(AT n )) 
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sup E 

nSN 



< oo. (2.23) 



Thus (2.12) holds. | 

(e) Proximity graphs. Devroye [7] defines a proximity graph on X to be 
one in which each {x,y} is included as an edge if a specified set S(x,y) is 
empty. If S(x, y) is the ball with opposite poles at x, y, then the associated 
proximity graph is the Gabriel graph. If S(x, y) is the intersection of B(x; \y— x\) 
and B(y; \y — x\) then it is the relative neighborhood graph. For a survey of 
applications of proximity graphs such as these, and also of the sphere of influence 
graph, see Jaromczyk and Toussaint [14]. 

Both the Gabriel graph and the relative neighborhood graph are translation 
and scale invariant, and stabilize on Vi, and also satisfy (2.12), subject to 
conditions on <f> and / similar to those already given for the sphere of influence 
graph; see remarks in [23], Section 9. Therefore Theorem 2.2 yields information 
about these graphs, adding to results in Devroye [7] on the expected number 
of edges. Numerical estimates for values of mean edge length and mean degree 
of various proximity graphs over homogeneous Poisson point sets are given by 
Smith (chapter III.C of [26]). 

(f) Boolean models. As already indicated, Theorem 2.1 extends to marked 
processes via the limit (2.8). An example of application is to Boolean models, 
whose importance in stochastic geometry and spatial statistics can be seen from 
e.g. Hall [10] and Molchanov [20]. Let us be a shape distribution, that is, a 
probability distribution on the space S of all compact sets in M. d . For measure- 
theoretic details see Matheron [17], page 27. Assume that /is is concentrated 
on connected sets contained in B(0;K) for some fixed finite K (i.e., uniformly 
bounded connected sets). On a suitable probability space let (Si, i > 1) be 
a family of random closed sets each with distribution fis, independent of each 
other and of (X\, X2, ■ ■ ■ ) (as usual, X\, X2, ■ ■ ■ are i.i.d. d- vectors with common 
density /). Let S n := U" =1 (A i + n _1 / d S' i ). We refer to Xi + n~ x / d Si as a random 
shape centered at Xi. The random set 5„ is a scale-changed Boolean model in 
the sense of Hall [10], pages 141 and 233. 

A connected component of S„ is often called a clump. A clump of order k 
is one which comprises precisely k random shapes. Let £/&(£„) be the number 
of clumps of order k and let [/(S„) := J^ fc £/fc(S n ), the total number of clumps. 
Let V"(S n ) denote the total volume of the set S n . In the case d = 2, we consider 
the total curvature functional of the set S„ (the product of 27r and the Euler 
characteristic of S„; see Hall ([10], Ch. 4.3)), which we shall denote W(lE, n ). 

Another statistic of interest is the off-line packing functional for the col- 
lection Xi + n~ x / d Si, 1 < i < n. This functional, denoted by M(S„), 
is the maximal number of non-intersecting random shapes in the collection 
Xi + n~ 1 ^ d Si, 1 < % < n. Additional statistics associated with Boolean models 
are found in Molchanov [20], for example. 

Let Sqo^ be the infinite Boolean model UxeVxiX + Sx), where V\ is a 
homogeneous Poisson process of intensity A on R d , each point carrying an in- 
dependent S- valued mark Sx with distribution /ig. Let A c := X c (d, us) denote 
the continuum percolation threshold, i.e., let A c be the supremum of the set of 
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values of A such that ^oo.a almost surely has no infinite connected component 
(note A c = oo if d = 1). 

Theorem 2.7 (i) There exist constants Ufc ]00 (k £ N) and Uoo (dependent 
on lis), such that n Uk{^ n ) — > Ufe ]OC in L 2 (for each k e N) and such that 
n~ 1 U(E n ) — » Uoo in L 2 . 

(ii) There exists a constant i>oo (dependent on /is) such that V{^ n ) — * ^oo 
in L 2 . 

(Hi) If d = 2 and the measure /j,s is concentrated on convex sets, then there 
exists a constant Woo (dependent on us) such that W(S n ) converges in L 2 to 

Woo. 

(iv) Ifsup xeK d f(x) < X c , then there is a constant m^ such that n Af(E n ) 
converges in L 2 to m^. 

Part (i) of this result adds to existing results in the literature such as Hall's 
result ([10], Theorem 4.7) on the number of clumps of order 1, especially since we 
do not restrict attention to a uniform underlying density / for the points X{. As 
usual, the value of the limiting constant u^ is of the form J E[£oo(V y ( x ))]f '{x)dx , 
where £(x; X) is described in the proof below. 

Part (ii) shows that the volume functional satisfies a weak LLN over non- 
uniform point sets, adding to results of Hall ([10] , Ch. 3.4) involving the vacancy 
functional of S„. Part (iii) is a weak LLN for the total curvature functional over 
non-uniform samples, and also adds to results of Hall ([10], Ch. 4.3)). 

In part (iii), the constant w^ is again of the form J E[£ :00 (Vf(^ x ))]f(x)dx, 
and if /Us is isotropic, there exist analytic formulae for ^[Coo('Pa)]; see (2.14) of 
[20], or (4.28) of [10]. 

Concerning part (iv), the off-line packing functional can be shown to be 
subadditive, and methods based on this fact show that a weak LLN also holds, 
at least for the uniform distribution on a cube of volume 1/A, even in the 
supercritical case A > A c . 

Sketch of proof of Theorem 2. 7. Define the rescaled Boolean model ~EJ n by 

Let the mark space be S with mark distribution fj,s- 

(i) Given a marked point set X C M. d with marks S x , x G X, let £(x; X) be 
the reciprocal of the order of the clump of U y< =x(y + S y ) containing x + S x . 
Then H^(X) is the total number of clumps of U y ^x(y + S y ); hence 

H in (X n ) = U(E' n ) = U(E n ). 

Stabilization of £ follows from the fact that fis is concentrated on uniformly 
bounded sets. The moments condition (2.5) follows from the uniform bound 
£(x;X) < 1. Therefore the LLN for U(E n ) follows from (2.8). The LLN for 
C/fe(S„) is proved similarly. 

(ii) This time let £(x; X) be the volume of the intersection of U y ^xS y with 
the Voronoi cell around x for X. Then 

H u (X n ) = U(E' n ) = nU(E n ). 
15 



Since us is concentrated on sets contained in -6(0; K), for any x £ X, the inter- 
section of U ye xS y with the Voronoi cell around x is contained in B(x; 2K), since 
any point lying outside B(x; 2K) but in U ye xS y must be closer to some point 
y £ X \ {x} than it is to x. Both the stabilization and the moments conditions 
in the marked point process version of Theorem 2.1 hold as a consequence of 
this, and (2.8) yields the LLN for V(E n ). 

(hi) This time let £(#; X) be the contribution of the random set S x to the 
total curvature of the union Uy^x^y This gives us H^(X n ) = W{EL n ). By 
convexity, £(x; X) is uniformly bounded, which gives us the moments condition 
(2-5). 

By the assumption of uniform boundedness of the random sets, if So is an 
independent random shape inserted at the origin, changes outside B(0; 2K) do 
not affect the contribution of Sq to the total curvature. Thus, for all t, £ is 
almost surely stabilizing on V T - This enables us to deduce the stabilization 
condition in the marked point process version of Theorem 2.1, and, we therefore 
deduce the LLN behavior. 

(iv) Let £(x; X) be either 1 or 0, depending upon whether or not x + S x is 
included in the maximal subset of non-intersecting shapes. If there are several 
such maximal subsets, choose one in an arbitrary deterministic but translation- 
invariant manner. Then H^(X) is the maximal number of disjoint shapes and 
H^ n (X n ) — M(E' n ). Since £ < 1, this gives the moments condition (2.5). 

If / < A c , then there is almost surely no infinite cluster in Vf( x ) 1 for we 
arc in the subcritical phase of continuum percolation. Since there is no infinite 
cluster, inserting a random shape at the origin thus almost surely changes the 
order of only finite clusters and thus the packing functional M stabilizes. | 

(g) Packing processes. Consider the following prototypical random sequen- 
tial packing model. Let B n ,i, B n ,2..., B n ^ n be a sequence of d-dimensional balls 
of volume n _1 whose centers are independent random d- vectors with common 
probability density function / : M. d — > [0, oo). Let the first ball B n ^\ be packed, 
and recursively for i — 2, 3, . . . , n, let the i-th ball B n ^ be packed iff B n ,i does 
not overlap any ball in B n _i, ...,B n ^-\ which has already been packed. If not 
packed, the i-th ball is discarded. Let Nf(n) be the number of packed balls, out 
of the first n to arrive. This is sometimes called 'on-line packing', in contrast 
with the off-line scheme described earlier. 

We may use our general result (2.8) for marked processes to obtain LLN for 
random sequential packing, as follows. 

For any finite point set X c M. d , assume the points have marks which are 
independent and uniformly distributed over [0,1]. Assume unit volume balls 
centered at the points of X arrive sequentially in an order determined by the 
associated marks, and assume as before that each ball is packed or discarded 
according to whether or not it overlaps a previously packed ball. Let £(x; X) be 
either 1 or depending on whether the ball centered at x is packed or discarded. 
Then with the binomial point process X n defined at (2.3), and £ n defined at (2.4), 
it can be seen that H^ n (X n ) has the same distribution as Nf(n). Following [24], 
we can show that £ is almost surely stabilizing on V T , t G (0, oo), with limit 
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£oo. Since £ is bounded it satisfies the moments condition (2.5) and therefore 
using the limit (2.8) we get the following LLN for Nf(n). 

Theorem 2.8 Let f : M. d — > [0, oo) be an arbitrary density. As n — > oo, 

n-'Nfin)^ f E[^ 00 {V f(x) )]f{x)dx in L 2 . (2.24) 

This result represents a finite input version of Theorem 5.1 of Penrose [22], but 
with a more general class of densities /. 

Theorem 2.8 extends to more general versions of the prototypical packing 
model. For example, by following the general stabilization analysis of [24], we 
can develop asymptotics in the finite input setting for the number of packed 
balls in the following general models: (i) models with balls replaced by particles 
of random size/shape/charge, (ii) time dependent, dynamic models, (iii) coop- 
erative sequential adsorption models, and (iv) ballistic deposition models. In 
each case, we obtain a LLN for the number of packed balls, from among the first 
n to arrive, when the distribution of the balls has a density / : M. d — ► [0, oo). 
See [24] for a discussion of these models and for laws of large numbers in the 
special case where the particless arrive uniformly at random over K d . 



(h) Combinatorial optimization. Both the shortest traveling salesman tour 
and the minimal matching are translation and scale invariant graphs on finite 
point sets. However, it is not known whether their definition can be extended 
to infinite sets in a manner that makes them stabilizing, and therefore we are 
at present unable to apply Theorem 2.1 or 2.2 in these cases. 

3 Proofs 

The proof of Theorem 2.1 centers around suitably coupling a version of the 
binomial process X n to a Cox process, that is, a Poisson process whose intensity 
measure is itself random. We do this as follows. On a suitable probability space 
suppose we have, independently, a d-dimensional variable X with density /, and 
a homogeneous Poisson processes V of rate 1 on M. d x [0, oo). 

Define coupled point processes V{n), A^_ 1 , and H n and a random variable 
£„, all in terms of V and X, as follows. Let V(n) be the image of the restriction 
of V to the set 

{(x, t)el d x[0,oo):t< nf(x)}, 

under the projection (x,t) i— » x. Then V(ri) is a Poisson process in R d with 
intensity function nf(-), consisting of N(n) random points with common density 
/. Discard (N(n) — (n — 1)) + of these points, chosen at random, and add 
[n— 1 — N(n)) + extra independent points with common density /. The resulting 
set of points is denoted X' n _ Yl and has the same distribution as <-f„-i defined 
earlier. 
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To define H n , let V n be the restriction of V to the set 

{(x,t):t<nf(X)}. 

Let Ti. n be the image of the point set V n under the mapping 

[x,t)^n 1/d (x-X). 

Given X = x, the point process V n is a homogeneous Poisson process of intensity 
1 on R d x [0,nf(x)], and therefore, given X = x, H n is a homogeneous Poisson 
process on M. d of intensity f(x). Define £„ to be the limit £oo (7i„). 

Then 7i„ is a Cox process, where the randomness of the intensity measure 
comes from the value of f(X). Note that the distribution of TL n , and hence that 
of £„, does not depend on n. 

Lemma 3.1 Given K > 0, we have 

lim P[n 1/d (X'_ 1 -X)n B{0: K) = H n n B(0; K)] = 1. (3.1) 

n — *oo 

Proof. Suppose X lies at a Lebesgue point of / (sec e.g. [25]). Given X = x, 
the expected number of points of V in B(x;Kn~ 1 / d ) x [0, oo) that contribute 
to H n but not to V{n) is 

n I {f{x) - f(y)) + dy 

J B(x;Kn- 1 / d ) 

which tends to zero because x is a Lebesgue point of /. The expected number 
of points of V in B(x; Kn~ 1 ! d ) x [0, oo) that contribute to V(n) but not to 7i n 
is 

n f (f(y) - f(x)) + dy 

JB{x;Kn- 1 / d ) 

which also tends to zero for the same reason. Finally the probability that 
V(n) n B(x;Kn^ l / d ) ^ X n _ x n B{x;Kn^ l / d ) tends to zero as n —> oo, since 
\N(n) — (n — 1)| is o(n) in probability. Integrating over possible values of X and 
using Dominated Convergence, we obtain (3.1). | 



Lemma 3.2 Suppose that £ is almost surely stabilizing on V T , with limit £oo (V T ), 
for all t e (0, oo), and £ satisfies the moments condition (2.5) for some p > 1. 
Then 



lim E[n^H u {X n )] = / E^ 00 {V f{x) )]f{x)dx. (3.2) 

n^oo J Rd 

Proof. Let e > 0. Then since 
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for any K £ N we have 

P[\UX;K-i)-Cn\>e] 
< P\n Y ^ d {X' n _ l -X)n B(0; K)^H n D B(0; K)] (3.3) 

+P[£(W„; K) - iiU n : K)>e\. (3.4) 

By the stabilization assumption, we can choose K > so that (3.4) is less than 
e/2, and then by Lemma 3.1, the expression (3.3) is also less than e/2 for n 
large. Since £„ have the same distribution for all n, it follows that 

UX;K-i)^Ci- (3-5) 

Since the bounded p-th moments condition (2.5) is assumed to hold for some 
p > 1, the variables ^ n (X; X^_ x ) arc uniformly integrable, and hence their 
expectations converge to that of (i. By conditioning on X, we obtain 

E[(i] = I E[U{V f(x) )]f{x)dx := p, (3.6) 

which is the right hand side of (3.2). Since E[H in (X n )] = nE[£ n (X; X^)], this 
gives us (3.2). | 

Next we consider E[£ n (Xi; X n )£ n (X2', X n )], and use a refinement of the cou- 
pling argument. Assume on a suitable probability space that we have, indepen- 
dently, two (i-dimensional variables X and Y with density /, and two homoge- 
neous Poisson proceses V, Q, both of unit intensity onl d x [0, oo). 

Define coupled point processes X' n _ 2 (a binomial process) and H.* , TC^ (both 
Cox processes) and variables C,f and (^f, all in terms of V, Q,X, and Y, as 
follows. Let X' n _ 2 be obtained just as X' n _ x was before, i.e. let V(n) be the 
image of the restriction of V to the set {(x,t) G M. d x [Ooo) : t < nf(x)}, 
under the projection (x, t) h^ x, and let N(n) be the number of points of V(n). 
Discard (N(n) — (n — 2)) + of the points of V(n), chosen at random, and add 
(n — 2 — N(n)) + extra independent points with common density /. The resulting 
set of points is denoted X' n _ 2 and has the same distribution as X n _2- 

Let Fx be the half-space of points in M. d closer to X than to Y, and let 
Fy be the half-space of points in R d closer to Y than to X. Construct Ti.„ as 
follows. Let V x be the restriction of V to the set F x x [0,nf(X)]; let Q x be 
the restriction of Q to the set Fy x [0,nf(X)]. Let 7i„ be the image of the 
point process V x U Q x under the mapping 

(x,t) >-> n 1/d (x - X). 

Given X = x, the point process V x U Q x is a homogeneous Poisson process 
of intensity 1 on R d x [0,n/(x)]. Hence, given X = x, Ti„ is a homogeneous 
Poisson process on M. d of intensity f(x); let C,f be the associated limit £oo(7~(-n)- 
Construct H.^ in the following analogous manner. Let V Y be the restriction 
of V to the set Fy x [0,nf(Y)]; let Q Y be the restriction of Q to the set 
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Fx x [0, nf(Y)]. Let H^ be the image of the point process V Y U Q Y under the 
mapping 

(x,t) *n 1/d (x-Y). 

By an argument similar to that used for W„ , the point process 7i„ , given Y = y, 
is a homogeneous Poisson process on M d of intensity f(y); we set (^ '■= £oc (TC^). 

We can now prove the following result, which is the case q = 2 of Theorem 
2.1. 

Proposition 3.1 If S, is almost surely stabilizing on V T , with limit £oo(^t); f or 
all t G (0, oo ) and if £ satisfies the moments condition (2.5) for some p > 2, 
i/ien as n ^ oo, 

n- l H^{X n )^ ( E[i oc {V f{x) )]f{x)dx ini 2 . (3.7) 

Proof. For any if > 0, we have 

hm P[n 1 / d («_ 2 U {F}) - X) n 5(0; X) = H* n B(0; if)] = 1; (3.8) 



lim P[n 1/d («_ 2 U {X}) - Y) n 5(0; X) = U\ n 5(0; 1H] = 1. (3.9) 

n — >oo 

The proof of these facts is just the same as that of (3.1) using the additional 
observation that 

lim P\B{X;KrT 1/d ) C F x ] = hm P[B(Y; &" 1/d ) C F Y ] = 1. 

n — >oo n — >oo 

Note that 7i% and Tt^ arc independent identically distributed Cox processes. 
Independence follows by conditioning on the values of X, Y; given these the point 
processes H.* and TC^ are constructed from Poisson processes on disjoint regions 
of space. Therefore, for each n, the variables (£ and (£ are independent. Also 
the joint distribution of C,^ , Cn is independent of n. 

Since by translation invariance 

UX;K-2 U {X,Y}) = mn l/d {{K-2 U {Y}) - X)), 
and 

UY; K-2 u {x, y» = e(o ; n 1 '^*^ u {x}) - r)), 

it follows from (3.8) and (3.9), by a similar argument to that which yielded 
(3.5), that asm oo, 

UX;K-2 u {x,Y})UY;<-2 u {x,r}) -^ cfCf, (3.io) 

and hence 

£„(Xi;#„)£ n (x 2; ;r n ) -H* cf Ci y - (3.n) 
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By assumption, the bounded p-th moments condition (2.5) holds for some 
p > 2. Then by Cauchy-Schwarz, 

supE[^ n (X 1 ;X n )^ n (X 2 ;X n ))P/ 2 } < oo 

and therefore the variables £ n (X\; X n )^ n (X2] X n ), defined for each n > 2, are 
uniformly intcgrablc, so that the convergence (3.11) also holds in the sense of 
convergence of means, i.e. 

lim E[i n {X 1 -X n % l {X 2 -X n )] = E[g&] = M 2 , (3.12) 

n — >oo 

with /i defined at (3.6) Here we have used independence of (£ and (^ . 
To complete the proof, observe that 

E^n^H^X^) 2 ] - n- x E\Z n {Xn X n ) 2 ] + (1 - i)£[£ n (X i; X n )^ n {X 2 ; X n )}, 

and in the right hand side the first term tends to zero by the bounded p-th 
moments condition, while the second term tends to /i 2 by (3.12). Therefore 
E[n- 1 H in (X n )] -*■ fi and £[(n _1 H in (X n )) 2 ] -> p 2 , so that n^H^jXn) tends 
to /i in mean square. | 

Finally we prove L 1 convergence, i.e. the case q = 1 of Theorem 2.1, com- 
pleting the proof of that result. 

Proposition 3.2 If '£ is almost surely stabilizing on V T , with limit £oo('Pt), f or 
all t G (0, oo) and if £ satisfies the moments condition (2.5) for some p > 1, 
then as n — > oo, 

n-'H^jX^^ f E[i oc {V f{x) )]f{x)dx ini 1 . (3.13) 

Proof. Given K > 0, define the functional 

( K (x,X) :=imn(^(x;X),K). 

Then by the stabilization condition for £, the truncated functional £ also 
stabilizes on V T with limit £%o(V T ) := niin^oo (V T ), K); we leave the reader to 
verify this assertion. Since £ K is uniformly bounded, by Proposition 3.1 and 
the fact that L 2 convergence implies L 1 convergence, we have 

n-'H^iXr,)^ [ E[£(V f(x) )]f(x)dx ini 1 . (3.14) 

JR d 



Moreover, 



< E^H^X) - rC 1 H tL K{X)\ 
E[^{n 1 l d X l -n 1 l d X n ) - S > K {n 1 l d X 1] n 1 ' d X n )] 
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which tends to zero as K — > oo, uniformly in n, because the assumed moments 
condition (2.5), p > 1, implies uniform integrability of the family of variables 
£(n 1 / d Xi;n 1/ ' d A' n ), n > 1. Also, by monotone convergence the right hand side 
of (3.14) converges to the right side of (3.13) as K — > oo. Combining these facts, 
and taking K to infinity in (3.14), we obtain (3.13).^ 

We now work towards a proof of Theorem 2.2. First we show that stabi- 
lization of the graph G as defined in Section 2.1 implies an apparently stronger 
form of the stabilizing property. 

Lemma 3.3 Suppose the graph G(X) is translation invariant and stabilizes on 
V\. Then with probability 1, for all points X G Vi.o there exists R(X) < oo 
such that the set of edges of G(Pifi) incident to X is unchanged if points are 
added and/or deleted outside B(0; R(X)). 

Proof. Let us say that a point iina locally finite set X is unstable for G(X) if 
there does not exist r > such that 

£{x; G{X)) = £{x; G(X n 5(0; r)) U A) 

for almost all finite A c K d \ B(0; r). 

Let t > 1. By translation- invariance and the fact that the Poisson point 
process is its own Palm distribution, the mean number of points of V\ in B(0; t) 
that are unstable for G(V\) is equal to 



/ P[x is unstable for V\ U {x}}dx 

JB(0;t) 



>B(0;t) 

which is zero by the assumption that G stabilizes on V\. Therefore, since t can 
be arbitrarily large, with probability 1 the Poisson process V\ has no unstable 
points. 

A further application of Palm theory for Poisson processes shows that the 
mean number of pairs of distinct points X, Y of V\ in B(0; t) such that Y is an 
unstable point of V\ , is equal to 

/ / P[y is unstable for V\ U {x,y}}dxdy 

JB(0;t) JB(0;t) 

and since the mean number of such pairs is zero, the above integral is zero, and 
therefore for almost all x in B(0;t), 

Leb {y g B(0; t) : P[y is unstable for V x U {x, y}}} > 0} = 0. (3.15) 

Choose xq in B(0; 1), such that (3.15) holds with x = xq. Then the set of z in 
B(—xo;t) such that 

P[z is unstable for V\ U {0, z}} > 

has zero measure. Integrating over z G B(—xo',i) and using Palm theory yet 
again, wc find that the mean number of Poisson points X G V\ in B(—xo',i) 
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such that X is unstable for G{V\ U {0}) is zero. Since t can be arbitrarily large, 
this gives us the result. | 

Proof of Theorem 2.2. To prove (2.10) let £(x; X) be the reciprocal of the 
order of the component containing x in G{X). If the component containing the 
origin in G(7 3 i,o) is Unite then by stabilization of G (Lemma 3.3) there exists 
R such that alterations to V\ outside B(0;R) will not cause any change in 
this component. If the component containing the origin in G(V\fi) is infinite 
then given e > we can find a connected subgraph of G(V\fl) of order greater 
than e _1 ; then by Lemma 3.3 there exists R such that alterations to V\ outside 
B(0; R) will not cause any removal of edges in this subgraph, and hence the 
functional £ stabilizes on V\ with limit a^ 1 . By scale invariance £ also stabilizes 
on V T with limit having the same mean as o~q . Since t;(x; X) is uniformly 
bounded by 1, the moments condition (2.5) is trivially satisfied and by Theorem 
2.1 we have (2.10). 

To prove (2.11) let £(x;X) be equal to 1 if G(X) contains a subgraph iso- 
morphic to r with a vertex at x, and to zero if not. This is bounded by 1, and 
by Lemma 3.3 £ stabilizes on V\ with limit equal to 1 if Eq occurs and equal 
to if not. Also £ stabilizes on V T with limit having the same distribution, and 
Theorem 2.1 applies to yield (2.11). 

To prove (2.13), observe that the functional L^(X) = YleeG(x) <?Kl e l) nas 
the representation L9(X) = L^(X) with 

i{x ] x) = 1 - J2 ^(N)- ( 3 - 16 ) 

ee£(x;G(X)) 

If G is scale invariant, then G(V T ) has the same distribution as G(r _1/ ' d 'Pi) and 
therefore £ stabilizes on V T with 

e\u{v t )] = \e J2 <H^ 1/d N)- ( 3 - 17 ) 

ee<£(0;G(-Pi,o)) 

Then (2.13) follows from Theorem 2.1. | 
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